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Abstract. We present the spline-interpolation approximate solution of 
the Dirichlet problem for the Laplace equation in the bodies of revolution, 
cones and cylinders. Our method is based on reduction of the 3D problem 
to the sequence of 2D Dirichlet problems. The main advantage of the 
spline-interpolation solution of the 3D Dirichlet problem is its continuity 
in the whole domain up to the boundary even for the case of the linear 
spline. 

1. Introduction 

We present the new — spline-interpolation — approximate solution of the 
Dirichlet problem for the Laplace equation which has wide applications to elec- 
tric and seismic prospecting [IJ [3J [H [51 [El H] • This method is applicable for 
the cylinders, cones and the bodies of revolution. There exist many methods 
of approximate solution of this problem. The classic approximate method is 
the finite-element one. The spline-interpolation method is the analog of finite- 
element method where the discrete boundary nodes are the closed boundary 
curves and the cells are the layers between two parallel planes. Our method is 
based on reduction of the 3D problem to the sequence of 2D Dirichlet problems. 
Note that one can apply the functions of complex variable in the case of 2D 
Dirichlet problem. There also exists another spline technique of solution, but 
the one presented here is substantially different from the one given, for exam- 
ple, in the book [2]. The main advantage of the spline- interpolation solution 
of the 3D Dirichlet problem is its continuity in the whole domain up to the 
boundary even for the case of the linear spline. The piece-wise analytic form 
of the spline-interpolation solution is also convenient for the applications. 

2. Formulation of the problem 

Let M be a body of revolution with the smooth generatrix x = <&(h) in 
XYH space. Let us denote by S the boundary of M . 

The classic internal Dirichlet problem for the Laplace equation 

(1) A 3 w = 0, 

is the following: given the function Uq on S it should be found the function 
u such that it satisfies equation (TTJ) everywhere in M and takes the given 
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boundary values Uq, that is 



u\s 



U . 



Here 




+ 




We call by the spline-interpolation solution of the problem the function which 
satisfies equation (TTJ) almost everywhere in M and takes the given continuous 
boundary values at a finite number of the curves being the sections of S by 
the planes h — hj. So we do not solve the initial boundary value problem but 
reduce it to the problem with the data given at some finite set of sections of S. 
For the case when the domain M has the upper and the lower ends Da and Dg 
on the planes h = A and h = B we assume that the given boundary values at 
these ends are polyharmonic functions. The spline-interpolation solution takes 
these given values at the ends. Note that the requirement of poly-harmonicity 
is not too restrictive because every continuous function given on the plane can 
be approximated with any accuracy and any smoothness by a polyharmonic 
function. 

It is well-known [8] that the external problem can be reduced to the internal 
one by Kelvin transform [7] . Let IR 3 \ M be the unbounded domain in which 
we must solve equation ([T]). We assume that the unit ball -Bi(O) is a subset 
of M. Then the inversion of R 3 \ M with respect to the unit sphere belongs 
to -Bi(O). It is well-known, that u(x, y, h) is the solution of equation (TTJ) inside 
the ball -Bi(O) if and only if the function v(x,y,h) — u(x/r 2 ,y/r 2 ,h/r 2 )/r, 
r = \J ' x 1 + y 2 + h 2 is the solution of ([U outside the ball B\ (0) which vanishes 
at infinity (O(i)). So we present and demonstrate the spline- interpolation 
solution of the internal problem for the body of revolution and only then discuss 
the case of the external problem and the problem for cylinders or cones. 



We construct the spline-interpolation solution for every fragment of M which 
is the layer between two planes h — hj, h = hj + \. The form of spline- 
interpolation solution is polynomial with respect to h. If we put the solution 
in the form 



in equation ([l) and equate the coefficients with the same powers of h we get 
the following relations: 



3. Analysis 



v 




(k + 2){k + 1H+2 + A 2 u k = 0, k = 0, . . . ,p - 2 



and 



A 2 itfc = 0, k = p,p- 1 



where 



^=d 2 x +d 2 y . 
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Now the solution u{x, y, h) satisfies equation ([T]) in the layer and u(x, y, h) is 
a ([§]+l)-haxmonic function for any fixed h. The main problem now is to satisfy 
the boundary condition and to glue continuously the solutions constructed 
for the adjacent layers. Such spline-interpolation solution is named here the 
continuous solution. The derivative Uh{x,y,h) of the continuous solution is 
bounded but discontinuous at the points of M on the cuts h — hj. We also 
show how to minimize these discontinuity. 

4. The continuous solution. The case of a domain without ends 

Assume that the body of revolution M lies between the planes h = A and 
h = B, the vertices Pa and Pb being the common points of the domain and the 
corresponding planes. Assume that the boundary values are given at m circles 
Cj, j — 1, m, which are the sections of S by the planes h = hj, j = 1, to, 
in the form fj : Cj — > R and also at the points Pa and Pb ■ 

We begin the construction of the spline-interpolation solution within the 
layer which contains one of the vertices, namely Pa, so it is the layer h G [A, hi]. 
Here we put 

v}-(x,y,h) = ul(x,y) + (h- hi)u\, 

where UQ(x,y) is a harmonic function in the domain which is the section of 
M by the plane h — hi and u\ is a constant. Evidently this spline satisfies 
equation (TT]). The function Uq{x, y) can be restored via the boundary values 
f\. We find the constant u\ using the relation 

i _ u(P A ) - u\{xa,Va) 
(A- ^) 

where u(Pa) is the given value at the vertex Pa, the parameters xa and yA 
are the abscissa and the ordinate of the point Pa, respectively. 
In the next layer h G [hi, /12] the spline has the form 

u 2 (x, y, h) = Uq(x, y) + (h- h^u^x, y), 

where both of the unknown functions Uq(x, y) and u\{x, y) are harmonic in the 
union of the projections of the sections of M by the planes h — h\, h = hi onto 
the plane XOY. Evidently this spline satisfies equation (fTJ. 

We put u^(x,y) = Uq(x, y). So the boundary conditions at the curve C\ are 
met. Now we find the boundary values of the function u\ which are equal to 

h - u 2 (C 2 ) 
h 2 — hi 

and restore the harmonic in the section of M by the plane h = h 2 function 
u\(x,y) via its boundary values. 
In the layer h G [J12, h^] we put 

u 3 (x, y, h) = ul(x, y) + (h- h 2 )u\(x, y), 

where the harmonic function u^{x,y) can be restored by setting Uq[x, y) = 
Uq(x, y) + (h 2 — hi)u\{x, y), and the harmonic function uf (x, y) can be restored 
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via the boundary values 

h - u 3 (C 3 ) 
h 3 - h 2 

We move from one layer to the next one successively and restore the functions 
Uq(x,u), Ui(x,y), k = 1, ...,m in the whole domain M. The construction for 
the last layer is similar to that for the first one. 

Note that the harmonic on one cut of M by h — hk function must be 
continuous and harmonic at the neighbor cuts. It holds true if the boundary 
data on the curves Cj, j — l,...,m, is given in the form of trigonometric 
polynomials. 

The constructed spline- interpolation solution is harmonic at the points of the 
common plane of two adjacent layers. The corresponding function is restored 
uniquely via its boundary values, so this solution is continuous over the axis 
OH. It is also continuous in any plane parallel to XOY up to the boundary of 
M. So the constructed linear with respect to h spline-interpolation solution is 
continuous in M up to the boundary. 

The derivatives of the constructed solution with respect to x and to y are 
continuous in M at every level h = const, but the first derivative with respect to 
h is discontinuous at the points of M at the levels h = hj, j = 1, m. We show 

that this discontinuity tends to zero when m — > oo and max \ hj+i — hj\ — > if 

j 

the boundary data of the corresponding classic Dirichlet problem is a smooth 
function and the boundary surface of M is smooth. 

Suppose that the boundary data of the classic Dirichlet problem is the 
function f(9,h), 9 € [0, 2ir], h e [A, B], which has the continuous deriva- 
tive fh(9, h), which belongs to Holder class H(a) with respect to h uniformly 
for 9 e [0, 2n]. So there exist the constants L > and a € (0,1] such that 
\f' h {9,h') - f h {9,h")\ < L\h' - h"\ a for any 9 e [0,2tt]. Let the generatrix 
x = < E > (/i) of the boundary surface of M be smooth: $'(/i) belongs to H(a). 
We assume that the distance between any two adjacent planes h = hj equals 
t = (B — A)/(m + 1). Consider two layers in M with the common plane h = 
and the splines constructed for these layers: 

u 3 (x,y,h) = u (x,y) +u{(x,y)h, h e [-t,0], 

and 

^ +1 (x,y,h) = u (x,y) +u{ +1 (x,y)h, h e [0,t]. 

We have 

u J '($(-t)cos6»,$(-t)sin6i,-t) = f(9,-t), 
m j '($(0) cos 9, $(0) sin 6>,0) = u J ' +1 ($(0) cos 9, $(0) sin 9, 0) = f(6,0), 
u j+1 {<f>{t) cos9,${t) sin9,t) = f(9,t), 

so 

u o ($(O)cos6>,$(O)sin0,O) = f(9,0), 
-tu{($(-t)cos9,$(-t)sm9) = [f(9,-t)-f(9,0)}+ 
+ [u ($(0) cos 9, $(0) sin 9, 0) - u Q {<f>{-t) cos 9, sin 9)], 

tu{ +1 {<$>{t) cos 9, <$>{t)sm 9) = [f(9,t)-f(9,0)} + 
+ [u ($(0) cos 9, $(0) sin 9, 0) - u ($(i) cos 9, sin 9)]. 
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We apply Lagrange formula and have 
u{(<f>(-t) cos 0, $(-*) sin0) = f t {6, ri{9)t)-[u' Qx {<5>{ii{9)t) cos 0, $(n(0)t) sin 0) cos 9+ 
W Oy (^(0)t) cos 6, $(iJ,(6)t) sin 6) sin 6}&(fi(6)t), 7?(0),M#) G (-1,0), 

cos0, sin 0) = / t '(0, £(0)t)-[Mo x ($Ou(0)i) cos 0, ®{n{8)t) sin9) cos6»+ 
+t^($(/i(fl)t)cosfl,$(i/(fl)t)sinfl)sinfl]$'(i/(e)t) ) e (0,1). 

Note that the functions u' 0x (x,y) and diffcrcntiable in the union of the 

disks with the radii $(—<), $(0) and 
Now 

uj($(0) cos0, $(0) sin0) = i*i ($(-*) cos6», $(-t) sin0)+ 

+[ui($(0) cos 0, $(0) sin 0) - uj ($(-t) cos 0, $(-t) sin 0)] , 

ui +1 ($(0) cos0, $(0) sin0) = cos 0, $(t) sin0) + 

+ K +1 ($(0) cos 0, $(0) sin 9) - u{ +1 cos 0, $(t) sin 0)]. 
Finally we have according to supposition 

K +1 ($(O)cos0,$(O)sin0) - uj($(0) cos0, $(0) sin0)| < 

The function y) — u\{x, y) is harmonic in the disk x 2 + y 2 < $(0), so 

\u{ +1 (x, y) — u{(x,y)\ < Kt a everywhere in this disk, and we can see that the 
discontinuity of the function Uh(x,y, h) vanishes when the distance t between 
any two adjacent planes h — hj tends to zero. 

We have the representation 

u{ +1 (<f>(0) cos 9, $(0) sin 9) - u{ ($(0) cos 0, $(0) sin 0) = t(j>{6, t) 

where <fi(9,t) is bounded when t tends to zero if the functions $"(ft) and 
fhh(0,h) are continuous. 

Note that we can construct the solution beginning from the lowest point as 
well as from the highest point. If U{(x, y, h) and Uf(x, y, h) are such solutions 
then aui(x, y, h) + (1 — a)uf(x, y, h) is also the solution for any a € (0, 1). 

5. The continuous solution. The case of a domain with the ends 

Here the domain M is situated between the planes h = A and h = B, the end 
Da is located on the plane h = A and the end Db is located on the plane h = B, 
the curves Cj are located on the planes h = hj, j = l,...,m,. The boundary 
data on the ends has the form UA{x,y), (x,y) G Ea and UB(x,y),(x,y) € Eb- 
The boundary values at the points of the curves Cj are given by the functions 
f(9,hj). 

The form of the solution depends on the order of harmonicity of the functions 
Ua(x, y) and Ub{x, y). Assume that the highest order of harmonicity of these 
functions is equal to (n + 1). 

We construct the spline for the layer h € [A, hi] in the form 

n 

u^x, y, h) = ul(x, y) + (h - A)u\(x, y) + u \k{x, y)(h - A) 2k , 

fe=i 
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where u\{x, y) — Ua(x,u). The spline satisfies equation (p} if and only if 
u 2/c(^ V) = A 2 «J(x, y), fc = 1, n, 

where A| is times applied Laplace operator A2 = d 2 + d 2 . Hence we find 
all coefficients with the even powers of h if we satisfy equation ([l} in the layer 
h £ [A, hi]. The only odd coefficient u\(x,y) must satisfy the equation 

A 2 u 1 1 (x,y) = 0. 

Therefore we have to restore the function u\(x,y) which is harmonic in the 
union of the projections of the sections of M by the planes h = A, h = hi onto 
XOY. We do it using the boundary condition 

For the next layer h £ [/ii, /12] we take the spline in the form 



u 2 (x, y, ft) = ul(x, y) + (h- hx)u\(x, y) + E uj k (x, y)(h - h 1 ) 2k , 

k=l 

where 

n 

u 2 (x,y) = ul(x,y) + ~ A)u\(x,y) + ^u\ k {x,y){h 1 - A) 2k , 

k—l 

and 

u\ k {x, y) = -j^-A k ul(x, y),k=l, n. 

The harmonic function u1(x,y) can be restored via the boundary data at 
the curve C%. 

All the splines for the next layers are constructed similarly. Only the last 
spline where we must satisfy both the given data on the section and on the end 
Db is different. 

We take the last spline for the layer h £ [h m , B] in the form 

u m (x, y, h) = vR(x, y) + (h- h m )u™(x, y) + 

n n 

+ E "Sfo y)( h - h ^) 2k + E u 2k+i( x > y)( h - M 2fc+ \ 



where 



fc=l k=l 



u% l (x,y) = u™ 1 (x,y) + (h m -h m _ 1 )u? 1 {x,y) + 



k=l 

t 
(2k)\ 



u" 2 l k (x,y)= { -^-A k 2 u^(x,y). 



So we glue continuously the last spline with the previous one at the points 
of the common section and provide equality of the coefficients with the even 
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powers of h in |T]). Now we must provide equality of the coefficients with the 
odd powers of h in |T]), that is why we put 

(-1)* 



< k+1 {x,y)= ^ k ^ iy A k 2 u?(x,y), k = l,..., 



Now equation ([1} is satisfied in the layer [h m ,B] and we must satisfy the 
boundary data at the points of Db- We write this boundary condition in the 
form 

- f-ll* 

E(^TT)I A ^ f(:K ' 2/)(jB " /lro)2fc+1 = 

= U B (x,y) - £ ^^(^KB " h m ) 2k+ \ 

where the right side of the relation is known (n + l)-harmonic function, we 
denote it by Fb(x,i/). Our aim is to restore the (n + f)-harmonic function 
u™(x,y) with the help of the identity 

" (—'\) k 

(2) £ -±—l—A k 2 u?(x,y)(B - M 2fe+1 = F B (x,y). 

k=0 { >' 

Let us apply the operator A2 to both sides of ([2j , now harmonicity of the 
function 

<+ifog)= (2n+l)l AS " r(a;>y) 

implies 

When we apply the operator A2 to both sides of the last relation, the func- 
tion u%! n _ 1 disappears from the left part. We continue application of the oper- 
ators till k — n. We obtain n additional equalities from equality After 
we apply the operator A% to <j2j) we have 

(In 4- 1 V 

and it is the expression of the harmonic function u™ n+1 (x, y). 
When we apply A^" 1 we have: 

(b - MArv + (b - /i m ) 3 Ar x < = 

= (B- / lm ) (2 _"~i ) 1 ! <_ 1 — (B — h m ) 3 ~A%u? = 

= (h m - B)^^^ -(B- hm )*^±±luZ +1 = Ar'Fs. 

After we put the known expression of u™ n+1 into the last relation we get the 
biharmonic function u^-i - We continue our inverse movement over the addi- 
tional equalities and restore all functions u™(x,y), beginning with u™ n+1 (x,y) 
up to u™{x,y). We get the last function directly from equality ©. When we 
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restore all functions u™ k+1 (x, y), k = 1, n, we have the last spline in the layer 

h e [/i m ,£]. 

The solution of the problem for a cylinder can be reduced to the solution 
for the circular cylinder due to the additional conform mapping. 

Note that as at the previous section we can construct the solution beginning 
from the lowest end as well as from the highest end. If ui(x, y, h) and Uh(x, y, h) 
are such solutions then aui(x, y, h) + (1 — a)uh(x, y, h) is also the solution for 
any a G (0, 1). 



6. Example 

Let the domain M be the half-ball bounded by the plane XOY from below 
and by the sphere x 2 + y 2 + h 2 = 1 from above. We construct the simplest 
spline-interpolation solution with the given data on the base h = 0, on the 
level h = 1/2 of the spherical surface and at the vertex at the level h = 1. We 
suppose that the data on the base is a biharmonic function on the variables x 
and y: 

N 

u(x, y, 0) = Re(J2(a° k - ib° k )(x + iy) k ) + 

fe=0 
N 

+ (x 2 +y 2 )Re(J2(al~ibl)(x + i y) k ). 

fe=0 



The data on the curve {x 2 + y 2 = |, h = j} has the form 

w(^cos0, ^-sm<t>, ^) = y + ^(flfe cos<i) + 6 fc sin(7)). 

k=0 



The data at the vertex is u(0, 0, 1) = u . 

We search for the solution for both layers h £ [0, |] and h e [|, 1] in the 
form of the polynomial of the second power over h. 
In the layer h G [0, 1/2] the spline has the form 

u(x,y,h) = ul(x,y) + hu\(x,y) + h 2 u\(x,y), 

where 

JV 

ul(x, y) = Re(^2(a° k - l b° k )(x + iy) k ) + 

N 

+^ 2 + y 2 )R<Y / (4-tb 1 k )(x + iy ) k ), 

k=0 
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N 

r a 



\{x, y) = 2[f + Re(J2(a k ib k ){-j=f{x + 



N 1 

-Re(J2(a k - l bl)(x + l y) k )-^ + 

k=0 ' 



1 N 

+ -Re(Y / (4-*b 1 k )(2k-l)( X + i y) k )}, 



k=l 
N 



ul(x, y) = -2[a\ + ite(£>J - ib\)(k + l)(x + iy) k )]. 

k=l 

In the layer h £ [1/2, 1] the spline has the form 

u(x, y, h) = ul(x, y) + (h- -)u 2 (x, y) + (h- -) 2 ul(x, y), 

where 

a N 2 

«<>(*, !/) = -£ + Rc(V(a fe - ib k ){^=) k {x + iy) k ) - 

JV TV 

--Re(J2(al ib\){x + zy) k ) + (x 2 + y 2 )Re(£{al - ib k ){x + ly ) k ), 

fc=0 fe=0 

5 

u\{x,y) = 2u + -aj - a , 



k=i 



N 

\k\ 



(x, y) = -2[a\ + Rc(J2(4 - ib\){k + l)(x + iy) h 



7. Approximation estimate of the continuous solution 

Assume that u is the difference between the exact solution of Dirichlet prob- 
lem and the spline-interpolation solution presented in the previous sections. 
This function vanishes at the points of Cj , at the points of the ends or at the 
top points. The function u(x,y,h) has continuous derivatives u' x (x,y,h) and 
u' y (x, y, h) at any intersection of the plane h — const with M. So we have the 
estimate 

\u{x,y,h)\ <u{h) + l max Ju£(x,y,h) + ug(x,y,h), 

where u(h) is the maximum of the module of the boundary difference between 
the classic solution and the constructed approximate solution on the level h = h, 
I is the distance from the point (x,y) to the boundary of the intersection of 
M with the plane h — h. Evidently, u(h) = for h = hj, j = 1, 2, m. The 
values of u(h) tend to zero when we increase the number of m and decrease the 
value max \hj+i — hj | if the corresponding classic Dirichlet problem has smooth 

3 

boundary values. 
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The estimate shows that the continuous solution is rather accurate in the 
neighbourhood of the boundary S of the body M. 

8. Construction of the smoothing solution 

Here we construct the continuous solution which minimizes the difference 

between the derivatives with respect to h of the adjacent splines at the common 

disk Ej. We consider a body of revolution without the ends or with one end. 

We suppose that the function is continuous and the boundary data of the 

l 

corresponding classic Dirichlet problem f(9, h) has the form Cfe(^) cos k9 + 

k=0 

Dk sin k9 where the functions Ck(t) and Dk(t) are at least 3-diffcrcntiablc, 
he [A, B\. 

The spline for the layer h <G [hj, hj+i] has the following form: 

2AT+1 

v? +1 {x,y,h) = u{ +1 {x,y)h\ 

where 

ui + p \x,y) = (-iyX>u^\x,y)/{2p)\, 

u{ + p 1 +1 {x,y)^{-lY^u{ +1 {x 1 y)/{2p+l)\, p = Q 1 ... 1 N. 

We choose the (7V+l)-harmonic coefficient ul +1 (x, y) so that u^ +1 (x, y, hj) — 
u?{x, y, hj), hence all even coefficients u 2p ~ 1 (x, y), p — 0, N, are known. We 
choose the (N + l)-harmonic coefficient u{ +1 (x,y) so that 

u j+1 (&(h j+1 ) cos6», &(h j+1 ) sm6, h j+1 ) = f(6, h j+1 ) 

and the norm of the difference (x,y,hj) — vP{x,y,hj) in the space L 2 is 
close to minimal. 

Let hj = 0, hj-i = -t, hj+i = t, &(hj) = 1, $(hj+i) = l + tb(t), where b(t) 
is bounded when t tends to zero. Gluing of the neighbor splines at the unit 
disk yields the equation u 2 p (x, y) = u 2p (x,y). 

After we consider the boundary values u^ +1 (x, y, t) and u^(x, y, —t) and ap- 
ply Lagrange formula we obtain in the same way as for the liner splines the 
representation 

N 

(3) EK^+i (cos0, sin 0) - ui p+1 (cos0, sin 9)}t 2 P +1 = t 2 ^(9,t), 

p=0 

where <j>(0, t) is bounded when t tends to zero. 
Consider the function 

JV 

u{x,y,h) = ^2[u2p+i{cos6, sinfl) — i4 P+1 (cosfl, sinfl)]/i 2p+1 = 

N 

= ^H2 P+ i(a ; ,y)/i 2p+1 

p=0 

in the intersection of M with the layer h € [0,i\. 
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Let us introduce the polar coordinates and take the coefficient u\(x,y) = 
u\ (r, 9) in the form of the following expansions 

N I 

ui (r, 0) = £ {af cos k6 + fi* sin k6)r k+2 P. 

p=0 k=0 

The last relation provides us with the expansions of the coefficients U2 P +i (r, 9) , 
p = 1, N. Due to representation ([3]) we can take the boundary data on the 
level h = t in the form 

I 

u($(t),0) = t 2 ^2(X k cos fc6> + Mfc sin fcfl) + t 3 ip(t,6). 

k=0 

Therefore we have the relation 
i i 
t[J2{a° k cos k9+ Pi smk9)(l+tb(t)) k +^2(a 2 k coskQ+ff. sink9)(l+tb(t)) k+2 +... 

fc=Q fc=0 

' f3 I 

+ J2( a l N cos ke + Pt N sin WO- + tb(t)) k+2N ] - ~E 4(A: + cos ^ + 

A:=0 ' fc=0 



A:=0 



+/3l smk9)(l+tb(t)) k + ...+Y^ 4N(k+N)(a 2 k N cos k9+p 2 k N sin k9)(l+tb(t)) k+2N - 2 } + 

t 2N+1 
(2N + 1)\ 



'' k=0 
I 

= t 2 ^2(X k cosk9 + [i k sink9) +t 3 ip(t,9). 

k=0 

The last relation yields the following representation of the coefficients a° k , j3% 
through the coefficients a? ,P k J ' , k — 0, .., I, j — 2, N, 

N 

(4) a° k = t\ k -a 2 -a*- ... - a 2N + t£4 a fc j + *V*(*)' 

3 = 1 
N 

(5) pl = tn k -pl-pl-...- p 2N + £ 4$ + t 2 Mt), 

3=1 

where <p k (t) and ip k (t) are bounded when t tends to zero. 

Now wc must choose the coefficients ctf?,^ so that the value ||«i||.l 2 (.e) 
becomes rather small. Here 

1 i 1 

IMlL(S) = 27T / [ag + air 2 + ... + a 2 N r 2N ] 2 rdr + ^[J K + 

k+1 

1 

+(&ir 2+k + ... + a 2N r 2N+k ) 2 rdr + [ {ft + p 2 r 2+k + ... + p 2N r 2N+k ) 2 rdr]. 
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Lemma 1 . The minimum of the quadratic form 
1 

[do + air 2 + ... + a N r 2N ] 2 rdr 

o 



N 

with the condition °>k = 1 is equal to 2 (n+i)' 1 • 

/c— 

Proof. We introduce the variable t = r 2 and represent the expression 
ao + ait + ... + a^t N as linear combination of polynomials Pkif) = Pk(2t — 1), 
t E [0, 2tt], where Pfc(t) are Legendre polynomials: 

N 



a + a x t + ... + a N t N = ^ b k Pk(t). 



This representation is unique and gives us the linear relations between the 
coefficients a k and b k , k = 0, ...,N. Now the given condition takes the form 

N N N 

fc=0 fe=0 k=0 

and the given quadratic form is equal to 

N 1,2 



1 P 1^2 

Q fc=0 fc=0 



We apply Lagrange method and find the values b k — ri|qjjp which provide 
minimum of the corresponding quadratic form with the given condition. So the 
minimum of quadratic form equals to 



IE 



1 % 



2^2fc + l 2(iV+l) 2 

We denote by a^ , k — 0, ...,N the values of the coefficients which provide 
minimum of the quadratic form given in the statement of Lemma 1. 

Now we put q.q P = tXotip , p = 1,...,N. So the first summand from © 
according to (Q| with k — has the following representation 

l l 

2ir J [a° + a 2 Q r 2 + ... + al N r 2N ] 2 rdr = 2vr j {t\ [l - (af + ... + ~a%) + 

o o 

+afr 2 + ... + d%r 2N ] + t 2 R (r, t, N)} 2 rdr = [1 + tRo(t, N)]. 

We need to find minima of some other quadratic forms but we have to 
introduce some polynomial bases different from Legendre polynomials in order 
to find the exact value of these minima. So we give here only the estimates of 
them. 
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Lemma 2. The minimum of the quadratic form 
1 



[do + r k { ai r 2 + ... + a N r 2N )] 2 rdr 



o 



N 

with the condition ak = 1 is less than ^3 

fe=0 

Proof. Wc have 



1 

a + r k ( ai r 2 + ... + a N r 2N )] 2 rdr = ^+2a J r k+1 (a ir 2 + ... + a N r 2N )dr+ 





1 

+ J r 2k+4 (a 1 + ... + a N r 2N - 2 ) 2 rdr. 




Let us put ao = 0. Now 
1 



J [a + r k ( ai r 2 + ... + a N r 2N )] 2 rdr = J r 2fc+4 (ai + ... + a N r 2N - 2 ) 2 rdr < 



< J (o a + ... + a N r 2N - 2 frdr. 


1 

The minimum of the quadratic form f(<Xi + ... + a,Nr 2N ~ 2 ) 2 rdr according 



to Lemma 1 is equal to 2^2- This quadratic form achieves the minimal value 
when ak = Oije—ii k = 1, ...,N. Lemma 2 is proved. □ 
Now we put 

Pi" = p = 1, N,k=l, I. 

and obtain according to (jlj and (0 the estimates of the following summands 
from ©: 

1 

- / («fc + «^ 2+fc + + «fr 2W+fc ) 2 rrfr < ^(1 + i? fc (t, iV)), 


1 

- / (# + /^ 2+fe + ... + /3fV 2Ar+fe ) 2 ^ < ^(1 + 0*(t, AO). 



Finally 



2 



where A is constant and $(t,N) is bounded for the fixed N when i tends to 
zero and therefore 

at 2 

(7) IK +1 -<IIL(^)<^ + ^(W 
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9. Approximation estimate for the smoothing solution 

In this section we present integral approximation estimate of the constructed 
solution. We suppose that the distance \hj + i — hj\ is unique for all j — 1, m, 
and equals to t = (B — A)/(m + 1). 

Assume that u is the difference between the exact solution of Dirichlet prob- 
lem and the smoothing spline-interpolation solution. Let us first consider the 
equality which holds true due to Gauss-Ostrogradsky theorem 

J J u(iii l dxdy+u x dydh+Uydhdx)+ S y^^ J J u(u : ' h (x,y,hj)—u : j^ 1 (x,y 7 hj))dxdy 

S 3 = 1 Ej 

= J J I l^ rac ^ :E ' ^' h)\ 2 dxdydh, 

M 

where u J (x,y,h) is the spline at the layer h € [hj-±,hj] and Ej is the disk 
which is the intersection of M with the plane h — hj. 
So 

III ^ rac ^' E ' ^' h)\ 2 dxdydh < max \u\ J J \gr&du\dS+ 

M S 
m p p 

+ max|u|^ / / \(u 3 h (x,y,hj) - u 3 h +1 {x,y, h 3 ))\dxdy. 

The right-hand side of the last inequality can be made arbitrary small due to 
the approximation of the boundary data and the smoothing of the derivatives 
with respect to h. Really 

Ef /■„ 1, . j+i, ... Dm(B-A) 

J J \(K( x >Vi h 3)- u h {x,y,hj))\dxdy < N ^ m + V} + 

(B — A) 2 - 
' (m + l) 2 

where D is constant and ^(t, N) is bounded when t tends to zero, due to ^7} 
and to Holder inequality. Let maxj;/ \u\ > 0. For arbitrary small e > we 
can choose N so large that the first summand in the right side of the previous 
inequality is less than e/(3max^/ |{t|). Now with the chosen number N we can 
choose the number m so large that the second summand in the right side of 
the previous inequality is less than e/(3maxM \u\). 

Obviously the value of maxs |u| [/ / |grad-u|(i5 can be made less than e/3 for 

s 

sufficiently large number m. 

Therefore J J J \gr&du(x,y, h)\ 2 dxdydh < e for large values of N and m. 

M 

Hence the Holder inequality implies 



5radu(x, y, h)\dxdydh < y / V(M)e, 
where V(M) is the volume of M. 



M 
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Now we prove the following estimate: 

J J \u(x, y, li)\dxdy < Q^fe 



for any h € [A, B], where E~ h is the section of M by the plane h — h. 

We prove this estimate at first for the case when M is a convex body of 
revolution. 

Let us fix h and choose a point s* G U, where U = Uj=i Q U Ca U Csi Cj> 
j = 1, n, are the curves with the given boundary data on S and Ca, Cb are 
the corresponding edges (or points). The cone CH(Ph,s*) which is the convex 
hull of Ph 1J{ S *} is a subset of M since M is convex. Let us consider the set 
{"f^x, y, s*)} of the straight lines, connecting the points (x,y) of P^ with the 
point s* and parametrized linearly by h: 

T h (x, y, s*) = (x(h),y(h), h), h g [h, h*], 

where h* is the /i-coordinate of s* . One can find the characteristic 

K(h,s*)= max 

Now we find the characteristic K{h) = min s * e ;j K(h, s*) and the point s*(h) 
which realises this minimum, the /i-coordinate of the point s*(h) being h**. 

The values of the function u vanish on the set U, so for any point p G P^ we 
have 

u( P ) = / d imf^ dh = J (grad ^, y> , ); 

h h 

Hence the estimate of implies 



i(p)| <#(/»)! J \\grsdu{x,y,h)\\dh\. 



Thus 



J J ' \u{x,y,h)\dxdy = J J J |(gradu(a;, y, h), ^j-)\dxdydh < 

P- h CH(P- h ,s*{h)) 

<K(h) J J J \gradu(x,y, h)\dxdydh < K(h)y/V(M)e = QVe. 



Here Q = V / V(M) max {K(h)}. 

he[A,B] 

For another types of bodies we can obtain similar estimate if we divide any 
section P^ into finite number of domains, construct the cone with the vertex 
from U which lies in M for every domain and take the integral of |u| over P^ 
as the sum of the integrals over the introduced domains. 
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10. On the similar problems 

As it was noted above the external problem for the body of revolution can be 
reduced to the internal one when the unbounded domain in which we solve the 
problem contains any ball in its exterior. We make the inversion with respect 
to the corresponding sphere and obtain the internal problem for the bounded 
domain. We construct the spline-interpolation solution of this internal problem 
and then make the inversion once more. 

The most simple case of the external problem is the case of the problem for 
the unit ball -Bi(O) exterior. The points of the unit sphere are invariant under 
the inversion. So the boundary conditions remain the same. 

We also can apply the spline-interpolation solution presented in this paper 
to the Dirichlet value problem for a cylinder or a cone with the section which 
is the known conformal map of the unit disk. Solution of 2D Dirichlet problem 
for the corresponding cuts can be reduced to solution of this problem in the 
unit disk. 
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